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showing asymptotic scaling. Indeed in the strong coupling region a quite large range 
of P values exists where the fundamental mass agrees, within about 5% on the square 
lattice and about 10% on the honeycomb lattice, with the continuum predictions in 
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I. INTRODUCTION 



Recent numerical studies of lattice two-dimensional SU(A^) x SU(A^) principal chiral 
models, with the standard nearest-neighbour interaction 

Sl = -2iV/3^ReTr[f/(x)[/t(a;+^)] , = 4^ ^ (1) 

have shown the existence of a scaling region, where continuum predictions for dimensionless 
ratios of physical quantities are substantially verified [0,0. The scaling region begins at 
relatively small values of the correlation length well within the expected region of conver- 
gence of strong-coupling expansion. Moreover by performing a variable change from the 
temperature T to 

Te = , = ]vk ' 

where E is the internal energy, one can find agreement in the whole scaling region between 
the measured mass scale and the asymptotic scaling prediction, within few per cent 0]. 

As a matter of fact, this may be considered as an evidence for asymptotic scaling within 
the strong-coupling regime, motivating a test of scaling and asymptotic scaling by strong 
coupling computations. As a byproduct, strong-coupling series can be analyzed to investigate 
the critical behavior of the N = oo theory, where Monte Carlo data seem to indicate the 
existence of a phase transition at finite (3. 

Ref. 1^ was devoted to a complete presentation of our strong coupling calculations per- 
formed by means of the character expansion. We calculated strong coupling series for several 
quantities on the square and honeycomb lattices. On the ordinary square lattice, we calcu- 
lated the free energy up to O and the fundamental Green's function 

= (lReTr[f/(x)f/(0)t]) (3) 

up to O For chiral models on the honeycomb lattice, defined by the nearest-neighbour 

action, longer series were obtained: the free energy up to O and G{x) up to O (/?^°). 
Lattice chiral models on square and honeycomb lattices are expected to belong to the same 
class of universality with respect to the continuum limit. As we will see from the strong 
coupling analysis, even at finite (3 large- chiral models on the honeycomb lattices show a 
pattern very similar to that observed on the square lattice. 

In this paper, which represents the logical continuation of Ref. [Q, we analyze the N = oo 
strong coupling series presented there and the results are compared with the continuum limit 
predictions and Monte Carlo simulations. The main result of our strong coupling analysis 
of 2-d N = oo chiral models on the square and honeycomb lattices is the identification of 
a scaling region where known continuum results are reproduced with good accuracy, and 
asymptotic predictions are substantially fulfilled in the energy scheme. 
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II. STRONG COUPLING EVIDENCE OF A LARGE-iV PHASE TRANSITION. 



Numerical simulations at large of SU {N) and U {N) lattice chiral models show evidence 
of a phase transition at iV = oo. Indeed sharper and sharper peaks in the specific heat 

^ = iVdT 

are observed with increasing N, suggesting a divergent large- limit at a finite /3 0. By 
extrapolating to iV = oo the positions of the specific heat peaks, we obtained a rather precise 
estimate of the critical coupling: Pc = 0.3057(3). Details on our Monte Carlo simulations 
and their analysis can be found in Ref. |^. 

In order to investigate the above issue, we analyze the N = oo strong coupling series 
by employing the integral approximant technique P,p!0|, which is especially recommended 
in the case of small critical exponent The method of integral approximants consists of 
representing the power series under study by the integral of a linear differential equation. 
In our analysis we considered the integral approximants obtained from a first order linear 
differential equation 

QMf'ix) + PK{x)f{x) + Rj{x) = 0(x^+^+^+^) , (5) 

where Ql,Pk and Rj are respectively L,K, and J order polynomials (we fix Qlq = 1). 
These approximants are singular at the zeroes Xq of Ql{x), and behave as 

A{x)\x- Xo\-^ + B{x) , (6) 

where A{x) and B{x) are regular in the neighbourhood of Xq, and 

^ _ PKixo) 

Given a M order series, L, K and J must satisfy the condition L + K + J + 2<M. 

Let us analyze the N = oo strong coupling series of the specific heat, which is even in /3: 

p-^C = 1 + 6 X + 30 + 266 + 2160 x^ + 19932 x^ 

+ 183638 x^ + 1754130 x^ + 16911192 x^ + 0(x^) , (8) 

where x = In Table | we report the first singularity in the real axis and the corresponding 
exponent for different values of L, K, and J. The results are quite stable, leading to a critical 
behavior of the specific heat typical of a second order phase transition: 

C ~ (9) 

From Table | we estimate 

/3e = 0.3058(3) , 

a = 0.23(3) . (10) 
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The errors are just indicative. They are the variance of the results in Table | after discarding 
the two furthest values from the corresponding average; they should give an idea of the 
spread of the results coming from different approximants. Notice that the strong coupling 
determination of /3c is in agreement with its estimate from numerical simulations at large 
N. 

As further check of the above resummation procedure, Fig. ^ compares, in the region 
(3 < (3c, SU{N) and U{N) Monte Carlo data of the specific heat at large N {N = 21, 30 
for SU{N), and = 15, 21 for U{N)) with the determinations coming from the resummed 
and the plain strong coupling series (|p. We recall that SU{N) and U{N) models should 
have the same large- limit. Monte Carlo data of C appear to approach, for growing 
A^, the determination from the resummed strong coupling series. As expected from simple 
considerations on the finite A^ corrections to the N = oo strong coupling series, U{N) models 
converge faster than SU{N) models to the N = oo limit in the strong coupling region. 

Monte Carlo data at large A^ seem to indicate that all physical quantities, such as the 
magnetic susceptibility x = J2x G{x) and the second moment mass {Mg = l/^G and 
X^Q = \ Yjx x^Gi^x) ), are well behaved functions of the internal energy even at A^ = oo 0. 
Therefore, as a consequence of the specific heat divergence, x and Mq should have a singular 
behavior with respect to [3. We indeed expect 



dlnx dlnM^ 
~d^ dp 



(11) 



in the neighbourhood of (3c- Notice that a behavior like Eq. (0) leads to a non-analytical 
zero of the /3-function Pl(T) at Pc- 

Pl{t) ^ \p-Pcr (12) 

around Pc, explaining the observed behavior with respect to P of the large- A^ Monte Carlo 
data for the fundamental mass 0. 

In order to check the behavior ([Til), we analyzed the corresponding strong coupling series 
by a modified integral approximant scheme forcing the approximant to have a singularity 



at /5 ^ 0.3058, obtaining biased estimates of the exponent in Eq. (|Tl|)- In this modified 
scheme the values of L, K, and J in Eq. (^) must be chosen according to the condition 
L + K + J+ l< M. We analyze the series 

=4 + 8p + 28p^ + ASP^ + 20Ap^ + AAOp^ + 17A0p^ + 3744/5^ + 15148/5^ + 35048/5^ 

dp 

+140980/3^° + 327600/3" + 1323612/3^^ ^ 3149112/313 ^ 12727908/3^^ + O (/3^^) , (13) 

and 

d In BM"^ 

tL^ = -4 - 10/3 - 28/3^ - 74/3^ - 224/3^ - 598/3^ - 1936/3^ - 5282/3^ - 17560/3^ 

dp 

-49170/3^ - 162144/3^° - 464426/3" - 1549656/3^^ - 4459234/3^^ + O (/3^^) . (14) 

In Table |I| we report the range of the exponent variations when varying the zero in the 
interval 0.3055 — 0.3061. These results are quite consistent with the exponent a obtained 
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in the analysis of the specific heat strong couphng series, supporting the relations ([TlD 
and therefore (p!2D . When performing an unbiased analysis of the series (13) and (|14|), 
that is without forcing the approximants to have a zero at a fixed f3, the singularity and 
the corresponding exponent turn out to be less stable; more terms in the series would be 
necessary to have a satisfactory analysis independent of that of the specific heat. 

In Fig. H we compare our strong coupling calculations of x with the Monte Carlo data 
of SU{N) and U{N) models at large A^. An improved strong coupling estimate of x? repre- 
sented by the full line in Fig. 0, was obtained by integrating the resummed series of d In x/ d/? 
(using L = 5, = 4 and J = 4, see Table 0). Growing iV, the SU{N) and U{N) data of 
X and approach the same N = oo limit, which is well reproduced by the resummation of 
the series (|13D and (|lj) in the strong coupling region. 

We should mention an apparent discrepancy between the SU{N) Monte Carlo results at 
large and the resummation of the N = oo series. In the SU{N) numerical simulations, 
when varying A^ (for sufficiently large A^) the position of the peak of the specific heat turns 
out to be quite stable with respect to the correlation length, ^^^'^ ~ 2.8, leading to the 
expectation that at the large- A^ critical point ^ 2.8 [0. On the other hand the estimate 
of from the resummation of the strong coupling series is slightly smaller: ~ 2.5. 



III. SCALING AND ASYMPTOTIC SCALING. 

In spite of the existence of a phase transition at A^ = oo, Monte Carlo data at large 
A^ showed scaling and approximate asymptotic scaling (in the energy scheme) even for 
P smaller then the peak of the specific heat [@]. The stability of this pattern suggests an 
effective decoupling of the modes responsible for the phase transition from those determining 
the physical continuum limit, and therefore that evidences of scaling and asymptotic scaling 
could be provided by the large- A^ strong coupling expansion. 

The on-shell fundamental mass M can be extracted from the long distance behavior of 
the correlation function in the fundamental channel G{x), or from the imaginary pole of 
its Fourier transform. We considered two estimators of M, fig and /id, defined from the 
long distance behavior of wall- wall correlation functions constructed with G{x) respectively 
along the sides and the diagonals of the lattice. An alternative mass Mg is defined from the 
inverse second moment of G{x). Unlike the on-shell mass M, Mg is an off-shell quantity, it 
is related to the zero momentum of the Fourier transform of G{x), indeed G{p)^^ ~ Mq+p^ 
at small momentum. 

The quantities /Xg, /^d and Mq enable us to perform tests of scaling based on rotational 
invariance at distances d ^ ^ = 1/ ^s, checking /is/ /id — 1, and on the stability of dimension- 
less physical quantities, looking at the ratio Hs/Mg- We should say that these tests concern 
the long distance physics of chiral models. 

Monte Carlo data at relatively large A^ showed that, within statistical errors of few per 
mille, the above scaling requirements are verified already at ~ — 2, well within the 
strong coupling region. Numerical simulations provided an estimate of the large- A^ limit of 
the ratio M/Mg in the continuum limit: M/Mg = 0.991(1) 0. 

In Ref. [|] the strong coupling series corresponding to the above mentioned quantities 
have been calculated, in particular Mq up to O (/?^^), = 2(cosh/is — 1) up to O (/9^^), 
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and Ml = 4(cosh/id/V2 - 1) up to O {13^^). 

In Fig. we plot the ratio yUg/ /^d vs. the correlation length = 1/-^g as obtained from 
our strong coupling series. The N = oo strong coupling curve confirms the large- Monte 
Carlo result: /is//^d — 1 within few per mille at ^ ~ 2. 

Fig. I shows the ratio fis/Mc vs. ^g- Notice the stability of the curve for a large region 
of values of C,g and the good agreement (well within 1%) with the continuum large- value 
extrapolated by Monte Carlo data. 

In order to test asymptotic scaling we perform the variable change indicated in Eq.(0), 
evaluating the energy from its strong-coupling series. The two loop renormalization group 
and a Bethe Ansatz evaluation of the mass A-parameter ratio ^ lead to the following large- 
asymptotic scaling prediction for the on-shell fundamental mass in the /Se scheme: 

M = ReAe,2iWe) , 
Re = 



Ae,2i{(^e) = ^J87tI3e exp{-87r /3e) , 

& = ^ . (15) 

In Fig. ^ the strong-coupling estimates of /is/A_B,2/ and Mg/Ae,21 are plotted vs. Pe, for 
a region of coupling corresponding to correlation lengths 1.5 ^ ~ 3. (We recall that 
Mg differs from M by about 1% in the continuum limit.) The agreement with the exact 
continuum prediction is within about 5% in the whole region. Notice also that both curves 
go smoothly through the value of Pe corresponding to the specific heat singularity /5c, which 
is f3^E^ ~ 0.220. 

The strong coupling curves in Fig. |^ were obtained from the plain series of the energy and 
respectively of and M^. In the case of Mg, we also determined Mg/A.e,21 evaluating the 
energy and by integrating the resummed series respectively of the specific heat and of 
d In /?M^/d/5. The resulting curve changes very little from that derived from the plain series, 
the difference between the two curves would not be visible in Fig. ^ This indicates once 
more that the change of variable P ^ Pe washes out the singularity in /3 when considering 
physical quantities. 



IV. CHIRAL MODELS ON THE HONEYCOMB LATTICE. 

On the honeycomb lattice we consider the action with nearest-neighbour interaction. It 
can be written as a sum over all links of the honeycomb lattice: 

= -2N(3 J2 Re Tr [Ui U}] , U G SU{N) , (16) 

links 

where /, r indicate the sites at the ends of each link. As on the square lattice, a lattice space 
a, which represents the lattice length unit, is defined to be the length of a link. The volume 
of an hexagon is fh = 3\/3/2. Straightforward calculations show that the correct continuum 
limit is obtained identifing 
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(17) 



A. The large- phase transition. 

On the honeycomb lattice we have calculated the N = oo strong coupling series of the 
specific heat, which is even in j3, up to 26*^ order in jS: 

p~^C = 1 + 10 + 90 + 396 + 728 + 9120 + 28186 x^ + 136800 x^ 

+886116 + 3129380 x" + +18935800 x^^ + q (x^^^ , (18) 

where x = (3"^. The integral approximant analysis of the above series, whose results are 



reported in Table III , leads again to a second order type critical behavior with the following 



estimates of the critical (3 and a exponent: 

(3, = 0.4339(1) , 

a = 0.17(1). (19) 

Notice that this estimate of the exponent a is very close to that of the square lattice. The 
uncertainty on both estimates cannot really exclude that they are equal, which would be an 
indication of universality. 

Also in this context we analyzed the strong coupling series of the logarithmic derivative 
of the magnetic susceptibility x ciiid (3Mq by the modified integral approximant method 
which forces the existence of a zero at (3c- In Table || we report the range of the exponent 
variations when varying the zero in the interval 0.4338 — 0.4340. As for the square lattice, 
the results in Table || are consistent with a divergence characterized by the specific heat 
exponent (cfr. Eq. (|TT|), supporting the existence of a non-analytical zero of the /3-function 
at Pc- 



B. ScaHng. 

On the hexagonal lattice the maximal violation of the full rotational symmetry occurs for 
directions differing by a 7r/6 angle, and therefore, taking into account its discrete rotational 
symmetry, also by a 7r/2 angle. So a good test of rotation invariance is provided by the 
ratio between masses extracted from the long distance behaviors of a couple of orthogonal 
wall- wall correlation functions constructed with G{x). 

In Ref . 0] we defined two orthogonal wall- wall correlation functions G^™^ (x) and G'^^ (x) , 
with the corresponding masses fii and fi2, which should both reproduce the on-shell fun- 
damental mass M in the continuum limit. In order to extract Hi and ^2 we evaluated the 
O (P^^) series of exp(— 3/ii/2) and the O (P^"^) series of exp(— \/3/i2/2). 

Fig. H shows the ratio yUi//i2 vs. = ^/Mq- As expected from the better rotational 
symmetry of the honeycomb lattice, rotation invariance is set earlier than for the square 
lattice: already at a correlation length — 0.5 /ii//i2 — 1 within 1%. 
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In Fig. 1^ we plot the ratio fii/Mc vs. ^g- The approach to the continuum hmit value 
seems to be substantially equivalent to that observed on the square lattice, but then for 

~ 1-5 the curve becomes unstable. Such an instability should be cured by an extension 
of the series. 



C. Asymptotic scaling. 

The asymptotic scaling test is again best performed in the energy scheme. This requires 
some weak coupling calculations, which present some subtleties on the honeycomb lattice. 
This is essentially due to the fact that, unlike square and triangular lattices, lattice sites are 
not characterized by a group of translations. Details on our weak coupling calculations are 
given in the Appendix. 

We calculated the internal energy (per link) up to two loops, finding 



N 6^3 



1 + — T= + o (t^) 



(20) 



The energy scheme consists in defining a new temperature Te proportional to the energy 

_ QV3N _ 1 

- N^l^ ' - iVT^ • ^^^^ 

The other important ingredient in this game is the mass A-parameter ratio in the hon- 
eycomb lattice regularization, which requires the calculation of the ratio between the A- 
parameter of the MS renormalization scheme and that of the honeycomb lattice reg- 
ularization Ah, given that the the (on-shell) mass A-parameter ratio in the MS scheme is 
known ||^. 

From a one loop calculation we obtained 

^ - Ae.J^^] (22) 

The ratio between Ah,^;, the A parameter of the energy scheme, and Ah is easily obtained 
from the two loop term of the internal energy: 

ej-^^]. (23) 



Then the N = oo asymptotic scaling prediction in the energy scheme is 



Rh,E = 8W — exp 



27r / 
— exn - 



TT 



e 1 3^3, 
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Fig. H shows the ratios ii\I^e,21 and MgI^e,21 vs. (5e (corresponding to correlation 
lengths 1 ^ ~ 2.5), as obtained from the corresponding strong couphng series. Again 
there is good agreement with the continuum prediction, especially in the region correspond- 
ing to correlation length ^ 2, where the agreement is within 10%. The curve correspond- 
ing to Mq is more stable, and it changes little when calculated resumming the involved 
series. 



APPENDIX A: WEAK COUPLING EXPANSION ON THE HONEYCOMB 

LATTICE. 

On the honeycomb lattice the sites cannot be associated to a group of translation. This 
causes a few subtleties in the analysis of models on such a lattice. 

The sites a; of a finite periodic hexagonal lattice can be represented in cartesian coordi- 
nates by 



h^x + km + k (^,o) 




(Al) 

We set a = 1, where the lattice space a is the length of a link. The total number of exagons 
on the lattice is L1L2, while the sites are 2L1L2. The coordinate k can be interpreted as 
the parity of the corresponding lattice site: sites with the same parity are connected by an 
even number of links. 

Notice that each of the two sublattices identified by x^{k,k) = x{k,k,—^) and 
x^{k,k) = x{k,k,^) forms a triagular lattice. Each link of the honeycomb lattice con- 
nects sites belonging to different sublattices. Triangular lattices have a more symmetric 
structure, in that their sites are characterized by a group of translations. It is then conve- 
nient to rewrite a field 0(x) = 4>{k,k,k) in terms of two new fields 0_(a;_) = and 
= 0(5+) defined respectively on the sublattices x_ and x+. A finite lattice Fourier 
transform can be consistently defined 



x± 

1 

VhLiL2 



x± 

^±{x±) = :;;rT-r T.^"''''^Mp) , (A2) 

p 



where = 3^/3/2 is the volume of an hexagon, and the set of momenta is 

^ 27r ^ 27r ^ 
p = —mipi + —m2P2 
Li L2 



mi = 1, ...Li , m2 = 1, ...L2 , 

'2 
.3 



Pi = (io) , p2 = (-^,^1 . (A3) 
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Notice that 

- - 27r 2-71 pi 

p-x = —hmi + — /2m2 + I3— ■ (A4) 
To begin with, let us discuss the simple Gaussian models, whose action can be written 

as 

= ^ E - <\>Mf (A5) 

links 

where xi^Xr indicate the sites at the ends of each link. Rewriting the field in terms 
of two fields 0_(a;_) and as described above, and performing the Fourier transform 

we obtain 

K 1 

-<p.{-p)<P+{p)H{-p) - , (A6) 

where 

H{p) = e-'P^^ (^1 + 2e*^cos^^ . (A7) 
From ( |A6|) we derive the propagators: 

a/3 1 

(0_(A;)0_(g)) = (0+(A;)0+(g)) = ^h— ^fc+g.o , 



where 



A(k) = - 
^ ^ 9 



2 — cos k2 cos —ki + cos / 



(A9) 



When x^ and x_ are the ends of the same link, i.e. |a;+ — a;_| = 1, one can easily prove that 

(M^^H-i^-)) - = ■ (Aio) 

The nearest-neighbour action of chiral models on the honeycomb lattice is 

= E 2ReTr [Ui t/t] , U G SU{N) , (All) 

links 

where /, r indicate the sites at the ends of each link. The perturbative expansion is performed 
by setting 
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u 



AA 



(A12) 



{Ta are the generators of the SU{N) group and Aa are A^^ — 1 real fields) and expanding U 
in powers of A. The action becomes 



^ j:\MAi-Arr + \Tr{Aj~AlY - ^TriA - Ar.)iAf - A^) + O {A<^) 



links 



(A13) 



The change of variables (|A12| ) requires the introduction of an additional term in the action 
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(AM) 



Then following the recipe illustrated in the Gaussian example, we rewrite the field A°'{xs) 
in terms of two new fields A°i{x^) and A^(a;+), whose propagators can easily derived from 
those of the Gaussian models, cfr. ( [ASI ). We are now ready to perform weak coupling 
calculations. 

Given the free energy per site 



■^Iny []d[/(x) exp(-Sh 



(A15) 



where Ug is the number of sites, the internal energy (per link) can be obtained by 



(A16) 



3 dp 

The internal energy up to two loops is given by Eq. (pOD. 

In order to evaluate the ratio between the A-parameters of the MS renormalization 
scheme and the honeycomb lattice regularization, we calculated the correlation function 



G{T,x^-y^) = -(ReTr[t/(x+)f/(i/+)t]) . 
In the x-space we obtained (neglecting 0(a) terms) 



GiT,x,a) = 1 + -^^TF{a/x) + O (t') 



where 

In the p space 

G{T,p,a) 



F{a/x) = ^ fin- -7s -In 2) 
27r V X / 



2N ^ 



1 + 



AN 



T[D{ap) 



3^3, 



O T 



(A17) 



(A18) 



(A19) 



(A20) 
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where 



D{ap) 



1 

2^ 



(lnap-21n2) . 



(A21) 



The above results required, beside the relation ( [AlOp , the calculation of the following inte- 
grals: 

dki dk2 e*^"^ — 1 



.2^ 27r J-^ 27r A{k) 

dki fit dk2 A{ p) - A{k) - A{k + p) 
A{k)A{k+p) 



F{a/x) + Oi- 



X 



2zL 27r J-^ 27r 

3 V3 



2D{ap) + 0{ap) , 



(A22) 



(A23) 



where the extremes of integration are chosen to cover the appropriate Brillouin zone, which 
can be determined from the finite lattice momenta (|A3|). 

The next step consists in determining the renormalized functions Z^'^^{T, afi) and 
Zif^{T,afi) that satisfy the equations 



T = Z^'^iT, a/i) t 



(A24) 



where t and Gj^ {t, x, fi) are respectively the coupling and the correlation function renor- 



malized in the MS scheme. In the MS renormalization scheme we have 



G^^^(t,a;/i = 2e~^^) = 1 + 0{t' 

t 



/i ' 2A^ p2 

Then by imposing Eqs. ( [A24D we obtain 



1 



Z^^{T,a^i) = 1 + T— (In a/i + rf) + O M 



where 



d 



-2 In 2 



2 27r 



iV2 3^3 ■ 

The constant d determines the ratio A;pg/Ah, indeed 

'N'^-2 2 



A, 



4exp 



TT 



(A25) 



(A26) 



(A27) 



(A28) 



For the interested reader we mention that Eqs. ( |A22[]A23| ) may be derived from the 
following exact result 



T" dki 



dko 



¥ 27r J-^ 27r A(A;) + m2 1 + 



— 1 + -m' 
2n\ 8 



-3/2 



1 + -m' 



-1/2 
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K 



1 H — m 
4 



1/2 



1 + -m' 



-3/2 



1 + -m' 



-1/2- 



(A29) 
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FIGURES 



FIG. 1. Specific heat vs. /5. The dashed and sohd lines represent the plain strong coupling 
series and its resummation. The estimate of the critical /? is indicated by vertical dotted lines. 
When error bars are not visible, they are smaller than the symbol size. 

FIG. 2. Magnetic susceptibility vs. (3. The dashed and solid lines represent the plain strong 
coupling series and its resummation respectively. The estimate of the critical f3 is indicated by the 
vertical dotted lines. 

FIG. 3. ^s/^d vs. = 1/Mg. 

FIG. 4. iXs/Mg vs. = 1/-^G- The dashed lines represents the continuum limit result from 
Monte Carlo data. 

FIG. 5. Asymptotic scaling test by using strong-coupling estimates. The dotted line represents 
the exact result ([l5|). 

FIG. 6. fJ'i/fJ-2 vs. = ^/Mg for the honeycomb lattice. 

FIG. 7. iii/Mg vs. ^g = for the honeycomb lattice. The dashed lines represents the 

continuum limit result from Monte Carlo data. 

FIG. 8. Asymptotic scaling test for the honeycomb lattice by using strong-coupling estimates. 
The dotted line represents the exact result (U). The full line corresponding to the ratio Mc/Ah^e 
was constructed by resumming the involved strong coupling series. 
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TABLES 



TABLE L Resummation of the strong coupling series of the specific heat. We analyze the series 
of /5~^C expressed in terms of cf. Eq. (^), for which M = 8. We report the first singularity in 
the real axis, /3o = \/xo, and the corresponding exponent versus L, K, and J. 
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TABLE XL Analysis of the series of dlnx/d/3 and dln/3M^/d/?. For some set of L,K and J 
we report the range of values of 7 corresponding to the range of zero values 0.3055 — 0.3061. 

L K J 7 

^ 4 5 4 0.20-0.24 

4 4 5 0.28-0.31 

5 4 4 0.23-0.27 

5 5 3 0.22-0.26 

6 3 4 0.22-0.26 
6 4 3 0.23-0.27 

4 4 4 0.23-0.26 

5 4 3 0.26-0.30 

5 3 4 0.07-0.12 

6 3 3 0.10-0.15 



dln(/3Mg.) 
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TABLE III. Resummation of the 24* order strong coupling series of the specific heat for the 
honeycomb lattice. We analyze the series of (3~'^C expressed in terms of cf. Eq. (|l8|), for which 
M = 12. We report the first singularity in the real axis, /3o = ^Jxq, and the corresponding exponent 
versus L,K, and J. 
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0.43495 
0.43312 


0.162 
0.165 
0.167 
0.171 
0.161 
0.163 
0.185 
0.171 
0.240 
0.101 


TABLE IV. Analysis of the series of dlnx/d/3 (19*'" order) and dln/3M^/d/3 (18*^ order) for 
the honeycomb lattice. For some set of L, K and J we report the range of values of 7 corresponding 
to the range of zero values 0.4338 — 0.4340. 
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